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>—> . ABSTRACT 

O ' 
t— i ; 

^ ' We consider the acceleration of energetic particles by a velocity shear in the 

t> ■ relativistic background flow containing scattering centers. Three possible acceler- 

^ ■ ation sites for astrophysical jets are identified: (1) gradual velocity shear parallel 

to the jet axis, such as a velocity profile decreasing linearly outward with radial 
' coordinates, (2) gradual velocity shear perpendicular to the jet axis, such as in- 

trinsic jet rotation, and (3) nongradual/ discontinuous, longitudinal velocity shear 
■ at the jet-side boundary. We determine the characteristic acceleration timescales, 

CLf specify the conditions for efficient acceleration, and discuss observational features 

with respect to each process. In particular, it is shown that in the case of (2) the 
higher energy emission is expected to be concentrated closer to the jet axis, while 
in the case of (1) and (3) the higher energy particles are likely to be located near 
the edges of the jet, thus possibly leading to some form of limb-brightening. 

C3 ' Subject headings: acceleration of particles — galaxies: jets 



1. Introduction 

Collimated relativistic outflows emerge from a variety of astrophysical sources ranging 
from active galactic nuclei (AGNs), //-Quasars, and neutron stars in binary systems, to 
gamma-ray bursts (Falcke et al. 1997; Zensus 1997; Mirabel & Rodriguez 1999; Fender et 
al. 2004). It seems that in all of these sources the relativistic outflows, or at least those 
parts observed in the radio band, are essentially launched from an accretion disk (Livio 
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1999; Marscher et al. 2002; Fender et al. 2004). In several sources, monitoring of their 
radio structures clearly indicates outflowing components at high (apparently) superluminal 
speeds (Mirabel & Rodriguez 1999; Jorstad et al. 2001), suggesting that the bulk plasma 
in these jets moves at relativistic speed along the jet axis, with the jets being observed at 
small viewing angles (Blandford et al. 1977). The acceleration of particles and the origin of 
emission in such jets has been widely analyzed under the assumption of a one-dimensional 
flow profile. While clearly quite instructive, such an approach appears adequate only to a 
first approximation, since real jets are naturally expected to exhibit a significant velocity 
shear across the jet. In the present paper we therefore aim to explore the relevance of 
different shear velocity profiles for the acceleration of energetic particles in relativistic jets. 
The paper is structured as follows. In § 2 basic results in the field of shear acceleration are 
reviewed. In §§ 3 - 5 three possible sites of particle acceleration in astrophysical jets (e.g., 
velocity shear parallel and transversal to the jet direction) are distinguished on the basis of 
theoretical and observational arguments. We analyze their associated acceleration potential 
using simplified shear flow profiles and determine their characteristic timescales. Finally, 
possible consequences are addressed in § 6. 



2. Physics and theory of shear acceleration 

Physically, particle acceleration in shear flows occurs as a consequence of energetic 
particles encountering different local flow velocities as they are elastically scattered off small- 
scale magnetic field inhomogeneities that are contained in the collisionless, systematically 
moving background flow. In each scattering event the particle momentum is considered 
to be randomized in direction but its magnitude is assumed to be conserved in the (local) 
commoving flow frame in which the electric field vanishes. Since the rest frame of the 
scattering centers is regarded to be essentially that of the background flow (i.e., second-order 
Fermi effects due to a random component of motion of the scattering centres are neglected), 
particles will gain neither energy nor momentum merely by virtue of the scattering if there 
is no shear (or rotation) present and the flow is not diverging. However, in the presence 
of a shear in the background flow, the particle momentum relative to the flow changes for 
a particle travelling across the shear. As the particle momentum in the local flow frame 
is preserved in the subsequent scattering event, a net increase in particle momentum may 
occur (Jokipii & Morfill 1990). 

The acceleration of energetic particles in collisionless plasmas with shear flows has 
attracted attention since the work of Berezhko & Krymskii (Berezhko & Krymskii 1981; 
Berezhko 1981, 1982). Starting with the nonrelativistic Boltzmann kinetic equation, they 
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showed that the particle distribution function in the steady state might follow a power law, 
/ oc p' ~^ +a \ with p' the magnitude of the commoving particle momentum, if the mean 
time interval between two scattering events increases with momentum according to a power 
law (r oc p' a , a > 0). The acceleration of cosmic rays at a gradual shear transition for 
non-relativistic flow velocities was studied independently by Earl, Jokipii & Morfill (1988). 
Assuming the distribution function to be almost isotropic in the commoving frame, they 
rederived Parker's equation (i.e., the transport equation including the well-known effects of 
convection, diffusion and adiabatic energy changes) but also augmented it with new terms 
describing the viscous momentum transfer and the effects of inertial drifts. Effects of an 
average magnetic field, not discussed by Earl, Jokipii & Morfill (1988), were included in 
the paper by Williams & Jokipii (1991), showing that the magnetic field could lead to an 
anisotropic viscosity if the particle mean free path is not sufficiently small (cf. also Williams 
et al. 1993). Jokipii & Morfill (1990) presented a microscopic analysis of the particle trans- 
port in a moving, scattering fluid which undergoes a non-relativistic, step-function velocity 
change in the direction normal to the flow. They demonstrated that particles may gain 
energy at a rate proportional to the square of the velocity change. Matching conditions in 
conjunction with Monte Carlo simulations for shear discontinuities were derived by Jokipii 
et al. (1989). 

Particle acceleration using Monte Carlo methods for a sharp tangential velocity discon- 
tinuity has also been studied by Ostrowski (1990, 1998). He showed that only relativistic 
flows are able to provide conditions for efficient acceleration, which could result in a very flat 
particle energy spectrum, depending only weakly on the scattering conditions. The possible 
relevance of such a scenario for the acceleration of particles at the transition layer between an 
AGN jet and its ambient medium was recently discussed by Ostrowski (2000) and Stawarz 
& Ostrowski (2002). Astrophysical applications of shear acceleration were also presented 
recently by Subramanian, Becker & Kazanas (1999); following earlier suggestions by Katz 
(1991), they investigated the bulk acceleration of protons in the tenuous corona above an 
accretion disk, assuming the magnetic scattering centers (i.e., kinks in the magnetic field 
lines with the lines themselves being anchored in the disk) to be driven by the shear of the 
underlying Keplerian disk. They demonstrated that shear acceleration may transfer the en- 
ergy required for powering the jet in blazar-type sources, accounting for bulk Lorentz factors 
up to ~ 10. 

The work on gradual shear acceleration by Earl, Jokipii & Morfill (1988) was success- 
fully extended to the relativistic regime by Webb (1989) (cf. also, Webb 1985). Assuming 
the diffusion approximation to be valid, he derived the general relativistic diffusive particle 
transport equation for both rotating and shearing flows. Subsequently, Green's formula for 
the relativistic diffusive particle transport equation was developed by Webb et al. (1994). 
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Applying their results to cosmic-ray transport in the galaxy, they found that the accelera- 
tion of cosmic rays beyond the knee by means of shear and centrifugal effects due to galactic 
rotation might be possible, but not to a sufficient extent. Using the relativistic transport 
theory developed by Webb (1989), Rieger & Mannheim (2002) devised a simple model for 
the acceleration of energetic particles in rotating and shearing relativistic astrophysical jets. 
By considering possible velocity profiles for such jets, they analyzed the effects of shear and 
centrifugal acceleration and showed that (1) in the steady state, power-law particle spectra 
may be formed under a wide range of conditions, and (2) the inclusion of centrifugal effects 
may lead to a flattening of the spectra. 

In this paper, we present for the first time an explicit derivation and analysis of the accel- 
eration timescales in relativistic gradual shear flows and quantify the conditions required 
for efficient particle acceleration. Observational and theoretical arguments are employed to 
distinguish between three possible acceleration sites and to show that shear acceleration of 
particles may be naturally expected in astrophysical jets. 

3. Longitudinal, gradual shear 

Generally, the interaction between a jet and the ambient medium is most likely to in- 
duce a velocity shear within the jet that is directed along the axis of propagation; of the 
type v = v z (r) e z . An energetic particle will experience a gradual shear profile if the ve- 
locity in the shear region changes continuously, and if its mean free path is much smaller 
than the transverse scale of the shear region. In the case of AGNs, for example, there is 
mounting observational evidence for jet stratification on the pc-scale, suggesting at least a 
two-component jet structure consisting of a fast inner spine (containing knots with trans- 
verse magnetic fields) and an extended, somewhat slower moving, shear boundary layer with 
a longitudinal magnetic field. Typical examples include the FR II source 1055+018 (At- 
tridge et al. 1999), the radio galaxies 3C353 (Swain et al. 1998) and M87 (Perlman et al. 
1999), and the gamma-ray blazar Mrk 501 (Edwards et al. 2000). The suggestions of a 
two-component jet structure are also supported by unification arguments showing that the 
radio and optical luminosities of BL Lac objects and FR I radio galaxies are not consis- 
tent with predictions from simple unification schemes, such as within a one-zone model. It 
was shown by Chiaberge et al. (2000), however, that the results can be reconciled within 
the unifying scenario if one assumes an internal velocity structure in the jet, where a fast 
spine (r spinc ~ 15) is surrounded by a slower, but still relativistically (Fi aycr ~ 2) moving, 
boundary layer. Relativistic hydrodynamical jet simulations may add further evidence re- 
vealing that the interaction between the jet and the ambient medium may indeed lead to 
a stratification of the jet in which a fast (T b ~ 7) spine or beam is surrounded by a slower 
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moving (r& ~ 1.7), high-energy shear layer (cf. Aloy et al. 2000; Gomez 2002). The shear 
layer in these simulations has been found to broaden with distance along the jet, reaching 
a transverse size comparable to the beam radius. There is also increasing evidence for jet 
stratification on larger scales. In order to obtain a good qualitative description of the total 
intensity and polarization systematics of the jets in a sample of FR I sources, Laing (1996) 
(cf. also Hardcastle et al. 1997) has suggested, that FR I jets consist of a fast spine with 
transverse magnetic fields, which decelerates from relativistic (3 S to nonrelativistic speeds on 
kpc scales, and a slower moving shear layer with a longitudinal or two-dimensional magnetic 
field, where the velocity decreases linearly with radial distance from the jet axis, which is 
probably caused by entrainment of external material. More specifically, it appears that the 
anticorrelation between brightness ratio (of jet to counter-jet) and width ratio from Gaussian 
fitting in the kpc-jets of FR I radio galaxies is suggestive of a relativistic flow profile that 
decreases radially outward from the jet axis, with the ratio /5 Zjmax / /3 Z;m i n of velocities between 
the center (the fast spine) and the slower moving shear layer (the edge of the jet) required 
to be in the range ~ (0.1 — 0.2) in order to produce the observed anticorrelation (Laing et 



In order to analyze the conditions for efficient particle acceleration in such shear flows, 
let us consider a simple background flow that moves at relativistic speed along the jet axis 
with magnitude depending on the radial coordinate r, so that in four- vector notation the 
bulk flow profile is of the form 



where p' is the commoving particle momentum and r ~ A/c is the mean time between 
two scattering events. For the above flow profile and by assuming isotropic diffusion, the 
relativistic shear coefficient T is given by (cf. eq. [A8]) 



al. 1999). 
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and hence the characteristic acceleration timescale in the commoving frame becomes 
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where v z = v z (r) and a power law energy dependence for the scattering timescale has been 
assumed, r oc p' a . 

Possible realizations of longitudinal shear flows include a velocity shear profile decreasing 
linearly outward (Laing et al. 1999) and a top-hat - type velocity profile as suggested from 
hydrodynamical simulations. The former can be modeled as a velocity profile with constant 
inner spine, i.e. relativistic r Zimax — c inside r\ and nonrelativistic v z ^ m at r 2 > r±, where 
the size of r 2 may be typically taken to be of the order of the jet radius so that 

Av 

v g (r) = v XtmBX --^(r-r 1 )H(r-r 1 ), (5) 

where H(r — r±) denotes the Heaviside step function, Av z = (v z>max — f z , m in) ~ c and 
Ar — {t2 — ri). The associated acceleration timescale is then of the order of 

3 (Ar) 2 
W(0 - -T TT^T ( 6 ) 



cA 7b (r) 2 l + 76 (r) 
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It seems clear that for efficient acceleration of electrons a shear layer of sufficient small radial 
extension is required. For typical relativistic, pc-scale radio properties in AGN jets [B ~ 0.01 
G, 7&(ri) ~ 10, 7 ~ 300 for radio emitting electrons and hence A ~ r g ~ 5 • 10 7 (7/3OO) cm], 
we may then have an electron acceleration timescale near 77 of order 

3 (Ar) 2 rnn /300\ / Ar V 
tsh,i ~ -^r^- ^ 500 yrs, 7 

% cr 9 V 7 / V0-003pcy 

while the corresponding isotropic electron synchrotron cooling timescale is of the order t coo \ ~ 
800 (3OO/7) (0.01 G/B) 2 yrs. Hence, for efficient electron acceleration one would need a Ar 
of the order of 0.001 pc or less. However, the conditions are much more favorable for proton 
acceleration, where (for producing emission at the same frequency, and all else being equal) 
Ar could be a factor m p /m e larger than the one required for efficient electron acceleration, 
so that for protons Ar may be of the order of 1 pc. This suggests that under general 
conditions shear acceleration of protons resulting in proton synchrotron emission is more 
likely to occur, even though usually the emission by protons may be neglected because of its 
relatively small emitted power. It could be shown that a similar conclusion holds for a top- 
hat velocity profile with a fast core region, an adjacent transition layer, and a slower outer 
collar extension, except that theoretically the transition layer no longer needs to be located 
close to the jet-side boundary. Generally, it is interesting to note that the shear acceleration 
timescale depends inversely on the particle mean free path, t sh l oc 1/A, or t shj i oc I/7 if we 
assume a gyroradius momentum dependence. The synchrotron loss timescale shows the same 
proportionality, i.e., t coo i oc I/7; hence once shear acceleration works efficiently, synchrotron 
losses are no longer able to stop the acceleration, and particles are continuously accelerated 
to higher energies unless a process like escape or cross-field diffusion becomes dominant. 
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4. Longitudinal, non-gradual shear 

We now consider a particle that becomes so energetic that its mean free path becomes 
larger than the transversal width of the velocity transition layer. The particle can then be 
regarded as passing unaffected through this layer and hence essentially experience a dis- 
continuous drop in velocity. If the particle elastically scatters off magnetic inhomogeneities 
frozen into the background flows and repeatedly crosses the layer, efficient particle accelera- 
tion to high energy may occur. Ostrowski (1990), for example, has argued that the interface 
between the jet interior in a powerful AGN and the ambient medium may represent such 
a transition layer across which a relativistic velocity drop naturally occurs. We therefore 
consider a velocity transition layer between two fluids A and B, where fluid A moves along 
the jet axis with velocity u\ and fluid B with velocity Relativistic particles with en- 
ergy Eq measured in the commoving frame A, which pass unaffected through the transition 
layer of width Ar, will have a total energy in the commoving frame B given by the Lorentz 
transformation 

E' = Ta Eq(1 + P cos 4>) ~ Ta Eq , (8) 

where cos0 = p- Au/[ \p\ ■ \Au\ ] with p the particle momentum, Au = (ui — U2)/(l — ui u 2 /c 2 ) 
the relative velocity difference, (3 = Au/c and Ta = 1/a/1 — (Au/c) 2 the associated Lorentz 
factor. The second relation on the rhs of equation (8) holds if the particle distribution is 
assumed to be almost isotropic near the velocity discontinuity. Hence, when a particle crosses 
the transition layer from A to B and elastically scatters off magnetic inhomogeneities frozen 
in the background flows, its energy has this value in the local rest frame B. If the particle 
recrosses the transition layer (from B to A), its total energy measured in the commoving 
frame A is then 

E = T A E'(l-{3cos ( j) , )~T A E' = T 2 A Eo, (9) 

where the second relation on the rhs again holds for an almost isotropic particle distribution. 
After two successive crossings, particles may thus increase their energy by a factor of T A . 
Hence, if the particle distribution is almost isotropic near the discontinuity, the mean energy 
gain for a single crossing of the discontinuity is given by (cf. eq. [8]) 

A E \ 

-g-^(r A -i). (io) 

It seems clear that according to this idealized picture, the increase in particle energy may 
be substantial provided the velocity difference is relativistic. For nonrelativistic velocities, 
however, only a mean energy gain of second order in Au is obtained. In the presence of 
relativistic flow speeds the full picture requires that the nonnegligible anisotropy of the 
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particle distribution is taken into account. The mean energy gain may then be expressed as 

/AM . . 

^J=Ve(T A -l), (11) 

with i] e < 1. Using Monte Carlo particle simulations within the strong scattering limit 
(AB/B ~ 1, kj_ ~ Ostrowski (1990) has shown that r\ e may still be a substantial 
fraction of unity. Technically, one can define a characteristic acceleration timescale given by 
^acc =< At > / < AE/E >, where < At > oc r is the mean time for boundary crossing and 
t = A/c, with A the particle mean free path. In the ambient medium rest frame one may 
thus finally arrive at (cf. Ostrowski 1998, 2000) 

4cc = a - c £ (1 - 10) ^ if r g > Ar , (12) 

where particle simulations suggest that a may be as small as ~ (1 — 10), provided that par- 
ticles are allowed to escape from the acceleration process once they have crossed a boundary 
at ?"max- Note that a is, of course, dependent on the radial escape boundary distance r max , 
so that for Au/c = 0.99 and r max /r s = 3, 9, 18, for example, one may have a ~ 4, 9, 18 (cf. 
Ostrowski 1990, Fig. 3a). In the presence of a longitudinal mean magnetic field, the particle 
mean free path may be of the order of the gyroradius r g (cf. Ostrowski 1998; Stawarz & 
Ostrowski 2002), but should be still larger than the width Ar of the transition layer, i.e., 
the acceleration mechanism may work fairly quick, provided we have very energetic particles 
with A ~ r g > Ar and r g < rj, where rj denotes the jet radius. Observations of pc-scale 
jets with evidence for a shear layer morphology (e.g., a boundary layer with parallel mag- 
netic fields) suggest that Ar < 0.5 r^. Taking Ar ~ 0.1 rj ~ 0.1 pc and B ~ 0.01 G for 
a semi-quantitative estimate, we would require very high seed electron Lorentz factors of 
about 7 e ~ 10 12 and proton Lorentz factors of about 7 P ~ 5 • 10 8 . The acceleration timescale 
should then be of the order of t acc ~ O.lpc/c ~ 0.3 yrs. Yet, due to the rapid radiation 
losses for electrons at the required high energies, the considered mechanism is unlikely to 
work efficiently for electrons since one typically expects r s (e) <^ Ar. Again the mechanism 
is much more favorable for protons; it may be quite possible to accelerate protons efficiently 
until their gyroradius becomes larger than the width of the jet. Note that for this case we 
have t acc oc A, while for the gradual shear cases above t s h,i oc 1/A. 



5. Transversal, gradual shear 

In addition to the possible occurrence of a longitudinal velocity shear in astrophysical 
jets, real jets are also likely to exhibit a velocity shear transversal to the jet axis. In par- 
ticular, several independent arguments suggest that jets may generally be characterized by 
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an additional rotational velocity component. The analysis of the jet energetics in extra- 
galactic radio sources for example has revealed a remarkably universal correlation between a 
disk luminosity indicator and the bulk kinetic power in the jet (Rawlings & Saunders 1991; 
Celotti & Fabian 1993) and established a close link between jet and disk (cf. also Livio 1999; 
Marscher et al. 2002), a fact which has been successfully exploited within the concept of a 
so-called jet-disk symbiosis (e.g., Falcke & Biermann 1995; Falcke et al. 1995). The results 
suggest that a significant amount of accretion energy, and hence rotational energy of the 
disk (cf. virial theorem), is channeled into the jet allowing for an efficient removal of angular 
momentum from the disk. Furthermore, direct evidence for intrinsic jet rotation has been 
recently established for stellar outflows (e.g., Bacciotti et al. 2000, 2002; Davis et al. 2000; 
Wiseman et al. 2001; Testi et al. 2002; Coffey et al. 2004). The optical jet in the T Tauri 
star DG Tau (central mass M* ~ 0.67 M Q ), for example, is associated with a Keplerian 
rotating circumstellar disk and appears to have an onionlike structure, with the faster and 
more collimated gas continuously bracketed by wider and slower material. The flow becomes 
gradually denser toward the central axis. High angular resolution observations revealed a 
systematic offset in the radial velocities of optical emission lines found in pairs of slits on 
alternate side of the jet axis, with average Doppler shifts being in the range of 5 — 20 km/s. 
This suggests a rotating jet flow with intrinsic, position-dependent, toroidal velocities in the 
range (6 — 15) km/s (for possible underestimation-effects, see however Pesenti et al. 2004) 
in the region probed by the observations, i.e. within several tens of AU from the jet axis 
(cf. Bacciotti et al. 2000, 2002; Testi et al. 2002). In the case of AGNs intrinsic jet rotation 
is also suggested by several indirect observations such as the helical motion of knots, the 
apparently oscillating ridge line in pc-scale jets, the detection of double helical patterns and 
the observations of periodic variabilities (e.g., Camenzind & Krockenberger 1992; Schramm 
et al. 1993; Steffen et al. 1995; Wagner et al. 1995; Krichbaum et al. 2000; Lobanov & Zensus 
2001). In the lighthouse model proposed by Camenzind & Krockenberger (1992), the origin 
of periodicity has been successfully related to differential Doppler boosting of emission from 
an off-axisymmetric density perturbation in the parsec-scale jet, which is dragged along with 
the underlying rotating bulk of the plasma. Finally, from a more theoretical point of view, 
intrinsic jet rotation is generally expected in magnetohydrodynamic (MHD) models for the 
formation and collimation of relativistic astrophysical jets (e.g., Begelman 1994; Camenzind 
1996; Fendt 1997; Ferrari 1998; Sauty et al. 2002). In such models, intrinsic rotation with 
speeds up to a considerable fraction of the velocity of light is a natural consequence of the 
assumption that the flow is centrifugally accelerated from the accretion disk. The set of 
available rotation profiles needs not necessarily to be disklike (i.e., Keplerian, cf. Blandford 
& Payne 1982), but may also include, for example, rigid and flat azimuthal velocity profiles 
(e.g., Vlahakis & Tsinganos 1998; Hanasz et al. 2000; Lery & Frank 2000). 
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In order to analyze the capability of particle acceleration by the shear in rotating flows, 
let us consider a simple model, in which the relativistic outflow velocity v z along the jet 
axis is kept constant and the velocity component perpendicular to the jet axis is purely 
azimuthal, described by the angular frequency profile Q(r). In the local commoving frame, 
the characteristic shear acceleration timescale for energetic particles with v' ~ c may then 
be estimated using (e.g. Webb 1989; Rieger & Mannheim 2002) 

tshear ^ = ^ ~ q) ^ f (13) 

where the shear coefficient T is given by (cf. Rieger & Mannheim 2002) 



1 / \4 2 

J^w) r 



dQ(r) 
dr 



1 2 



:i-v 2 z /c"). (i4) 



Here 7&(r) = (1 — f 2 /c 2 — f2(r) 2 r 2 /c 2 ) -1 / 2 denotes the position-dependent Lorentz factor of 
the background flow, A = A(7) = Ao7 Q < r g , with a > 0, is the particle mean free path, and 
7 is the Lorentz factor of the particle. Since the jet flow is likely to be launched from an 
accretion disk, one expects the jet flow to be radially confined to a region between rj n and 
r out ~ Ty Keplerian or flat rotation belong to the most relevant realizations of azimuthal 
shear profiles, and we will thus center our discussion here on these two types: 
(i) For a Keplerian rotation profile of the form Q(r) = Cl k (r in /r) 3//2 , where Cl k is a constant 
and r in < r < rj, Eq. (13) results in 

t k 20 c 1 1 If r \ 3 

shearW ^ 3 (4 + ft) A % 4 (l _ y 2 / ^ ^ [ r J " ^ > 

For a quantitative estimate we assume that because of a net radial mass flow inward, the 
rotational profile remains similar when the jet boundaries open slightly with distance, so 
that the rotational velocity at the inner boundary of the jet is of the order of the rotation 
velocity near the innermost stable orbit r is of the underlying accretion disk. We then have 
Clk — Vk(r[ n )/r is with r is = kor s , where r s denotes the Schwarzschild radius, and ko = 0.5 in 
the Kerr BH case, so that the acceleration timescale becomes 

k 20 c 1 1 1 2 2 / r \ 3 

WW ^ g-^^y j ~, (1 _ v2Jc2) k r s \ -J . (16) 

Using physically reasonable parameters, i.e., r in ~ 0.01 pc, v k (j-- m ) — £ k c with £ k < (1 — 
^/c 2 ) 1 / 2 ~ 0.3 for v z /c ~ 0.95, k = 0.5 (Kerr BH), and A(t) =t)r 9 ,r 9 ~ 5 x 10 7 (7/300) 
cm, 77 < 1, a — 1, r s ~ 3 • 10 13 cm (i.e. M B h ~ 10 s M ) and 7^ ~ 4, one can obtain an order 
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of magnitude estimate for the acceleration of radio emitting electrons in pc-scale jets given 
by 

4e» - 30 (™ )(-^— V yrs. (17) 



7 J \0.01pc 

The corresponding (isotropic) electron synchrotron cooling timescale, on the other hand, is 
of the order t coo i ~ 800(300/7) (0.01 G/B) 2 yrs. It seems clear that (a.) for electrons the 
acceleration process is rather limited and confined to a small inner region of the jet, while 
it is again much more favorable for protons, where (for producing emission at the same 
frequency, and all else being equal) tg hear p (r) ~ 0.7 (r/0.01 pc) 3 yrs, while t C ooi, P ~ 3 • 10 7 yrs, 
and (b.) substantial rotation is required for efficient acceleration. 

(ii) For a flat rotation profile of the form f2(r) = f2/ (r in /r) on the other hand, where f2/ is 
a constant and r in < r < rj, equation. (13) becomes 

/ 15 c 1 1 1 / r \ 2 

WW ^ j ~, (l-v*J<?) W f \~J ■ ( } 

For a quantitative estimate, we again suppose maintenance of the rotational profile, so that 



Vlf ~ v^j/ris where v^j = £/c, with £f < y 1 — v 2 /c 2 ~ 0.3. Using the same parameters 
as above (i.e., r; s = kor s with ko = 0.5, 7/ ~ 4 etc.), we may obtain an order of magnitude 
estimate for the acceleration of electrons by shear given by 

The acceleration potential is then somewhat improved compared to the Keplerian case, but 
it is still rather limited for electrons, while it is again much more favorable for protons. 
Our present analysis is aimed toward a study of the shear acceleration characteristics, and 
hence we have neglected centrifugal acceleration effects in rotating flows. We note however, 
that such effects may become particularly relevant for flat rotation profiles, in which the 
centrifugal acceleration timescale could be of the same order of magnitude as the shear 
timescale (cf. Rieger & Mannheim 2002). Generally, the shear acceleration timescale for 
rotating flows increases with the radial coordinate (i.e., with the cube for the keplerian case 
and the square for the flat case), and hence the higher energy emission is naturally expected 
to be concentrated closer toward the jet axis. As the synchrotron loss timescale scales with 
I/7, synchrotron losses are again no longer able to stop the acceleration once the shear 
mechanism has started to work efficiently (cf. § 3). 



6. Discussion and conclusion 



In the present paper we have identified three possible velocity shear sites in powerful as- 
trophysical jets on the basis of theoretical and observational arguments. It seems intuitively 
clear that the scattering of energetic particles across such shear profiles should naturally 
lead to a net increase in particle momentum, i.e., to the acceleration of particles. Using 
an idealized approach, we have specified order-of-magnitude estimates for the characteristic 
acceleration timescales using reasonable parameters for AGN jets. A more detailed determi- 
nation of the acceleration efficiency, which seems clearly desirable, however would require a 
detailed, but currently inaccessible knowledge of the real intrinsic jet properties and a more 
rigorous treatment of the physical mechanism responsible for the scattering. The last point 
is particularly illustrated by the dependence of the acceleration timescale on the particle 
mean free path (inversely proportional for the gradual shear cases and proportional for the 
nongradual shear) and, in the presence of a longitudinal mean magnetic field, by the need 
of a sufficiently strong diffusion perpendicular to the field. Our present approach for the 
gradual shear, for example, is based on an idealized large angle (hard sphere) scattering 
model, where n±_ ~ k\\ in the limit WjT < 1. While we thus caution that more detailed 
physical models need to be developed before further conclusions can be drawn, we never- 
theless suggest that the generic nature of our analysis might give valuable insights into the 
acceleration of particles by the shear in astrophysical jets. 

In general, the particle mean free path for electrons is much smaller than the particle mean 
free path for protons emitting at the same frequency. Shear acceleration is thus much more 
favorable for protons than for electrons. For the considered realizations, balancing the accel- 
eration timescales with synchrotron losses indeed suggests that shear acceleration is rather 
marginally significant for electrons, while it could be reasonably efficient for protons. 
Phenomenologically, one expects that particle acceleration by a (nongradual) shear at the 
jet-side boundary and electron acceleration by a (gradual) longitudinal velocity shear within 
the jet, which decreases linearly with the radial coordinate, are associated with the higher 
energy particles being concentrated closer toward the edges of the jet, and thus possibly 
related to some form of limb brightening given suitable flow velocities, jet orientations, and 
magnetic field strengths. The picture is, however, generally more complex if gradual shear 
acceleration of protons or top-hat velocity profiles are considered. For example, if a top-hat- 
type velocity profile with a radially extended collar is indeed physically realized, as perhaps 
suggested by two-flow models (Sol et al. 1989; Roland et al. 1994), the situation may be 
reversed, such that higher energy particles are located closer toward the axis. On the other 
hand, in the case of an intrinsically rotating flow, the shear acceleration timescale usually 
increases strongly with radial coordinate, so that higher energy particles are generally ex- 
pected to be concentrated closer towards the jet axis. It is, however, again possible, that the 
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real, observed situation is more ambiguous. In the case of an internally rotating jet flow, for 
example, a transition layer from the rotating jet interior to the nonrotating environment is 
still necessary. For a comparatively weak magnetic field and a sufficiently strong rotationally 
sheared layer this may lead to the onset of a strong magnetorotational instability (cf. Balbus 
& Hawley 1991), resulting in the growth of small-scale, tangential magnetic field discontinu- 
ities until the magnetic energy is eventually dissipated through reconnection events in which 
efficient particle acceleration may occur. In that case we may also have a contribution of 
highly energetic particles from the edges of the jet. 

So far, we have not considered the relevance of second-order Fermi effects that may con- 
tribute to the acceleration of particles (for a complementary approach, see Manolakou et 
al. 1999; Stawarz & Ostrowski 2002). Yet, if the scattering centers are not assumed to be 
completely frozen in the background flow, second-order Fermi effects should become domi- 
nant at lower energies. Quantifying these effects by determining the intrinsic Alfven velocity, 
for example, is, however, notoriously difficult because of our ignorance of the real physical 
parameters in astrophysical jets. Nevertheless, it seems qualitatively evident that, for a gyro- 
radius momentum dependence of the particle mean free path, the timescale for second-order 
Fermi acceleration scales as t2ndF oc A oc 7. The acceleration timescales in the gradual shear 
cases on the other hand, scales as t s h oc 1/A oc I/7. Hence, even if the observed lower en- 
ergy emission is mainly due to particles accelerated via a second-order Fermi process, there 
exists an energy scale above which gradual shear acceleration becomes dominant. In that 
case, second-order Fermi acceleration may be considered as providing the seed particles for 
efficient gradual shear acceleration within the jet. Accordingly, one may perhaps develop a 
picture in which supra-thermal particles, injected and centrifugally accelerated at the base of 
a jet (e.g., Rieger & Mannheim 2000), are - in the absence of strong shocks - initially further 
accelerated along the jet via second-order Fermi and then subsequently via gradual shear 
effects, until eventually (nongradual) shear acceleration at the jet-side boundary becomes 
possible, potentially allowing for the formation of a cosmic-ray cocoon (cf. Ostrowski 2000). 

We are grateful to M. Ostrowski and K. Mannheim for discussions and comments. 
Useful comments by an anonymous referee, which helped to improve the presentation, are 
also gratefully acknowledged. FMR acknowledges support through a Marie-Curie Individual 
Fellowship (MCIF-2002-00842). 
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A. Longitudinal, gradual shear flow 

In the commoving frame, the acceleration coefficient due to a velocity shear can be cast 
into the form (e.g., Webb 1989, eq. 3.27) 

AP ' V 1 V'tF) (Al) 



At / p' 2 dp' 

where p ' denotes the commoving particle momentum, with p ' ~ p ' for the energetic parti- 
cles considered here, r ~ A/c is the mean scattering time,- and T is the shear coefficient. In 
the strong scattering limit (i.e., w 9 'r< 1, with oj g ' the relativistic gyrofrequency measured 
in the commoving frame), we can use (see Webb 1989, eq. 3.34) 

r = ^°°e° aP , (A2) 

where a a p, with a, (3 — 0, 1, 2, 3, is the (covariant) fluid shear tensor given by 

2 c 
e a p = V a u p + V p u a + u a U[3 + iipu a + - (g af3 + u a up) V ' &u . (A3) 

Here g a p denotes the (covariant) metric tensor, e.g., for cylindrical coordinates x a = (ct, r, <p, z) 
we have (g a /3) — diag(— 1, l,r 2 , 1) (e.g., Rieger & Mannheim 2002), and V Q denotes the co- 
variant derivative. Note that for cylindrical coordinates, the only non-vanishing connection 
coefficients (Christoffel symbols of 2nd order) are 

Tl 2 = -r, 1^ = 1^ = 1. (A4) 
For a time- independent, relativistic bulk flow velocity profile of the form 

u a = lb (l, OA—) , (A5) 



with 7b(r) = 1/a/1 — v z (r) 2 /c 2 the bulk Lorentz factor and u a = 75 (— 1, 0, 0, v z (r)/c) both, 
the fluid four divergence, i.e., Vpu^, and the fluid four acceleration ii a = u^V^Ua vanish. 
Moreover, it is easy to show that the only non-vanishing components of the shear tensor are 
given by 

a„ = = =g£- (A6) 
a 13 = ^ = ViU 3 = — = -h b (r)— (A7) 
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Noting that a i = — o" 01 and a 13 = a 13 the relativistic shear coefficient could then finally be 
written as 



dv z (r) 
dr 



v z (r-y 



(A8) 



which, for nonrelativistic flow speed (i.e., 7^ — > 1 and v z /c — > 0) reduces to the nonrelativistic 
shear flow coefficient (e.g., Earl, Jokipii & Morfill 1988). For an energy-dependent scattering 
timescale of the form r oc p ' a the characteristic acceleration timescale can be expressed as 



W( r ) 



p 



15 c 



1 



<P>> (4 + a) A 7fc(r)2 (^r 1 + 7k(r)2 ^ 



(A9) 



where i> 2 = i> 2 (r), and where in the presence of a background magnetic field along the jet axis 
A = A(7) has to be smaller than the gyroradius to allow for isotropic diffusion characterized 
by k± ~ km, where (e.g., Webb 1989, eq. 3.7) 



K\\ 



1 , 

- V T c 

3 



and 



K\ 



K± 



;i+^ 2 r c 2 ) 



(A10) 



denote the parallel and perpendicular diffusion coefficients. 
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